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Prediction  of  annual  sunspot  numbers 
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*  1«  Purpose 

Since  discovery  of  the  sunspot  cycle  by  Schwabe  in  1851,  conside¬ 
rable  interest  has  been  exhibited  by  a  number  of  investigators  in  fore¬ 
casting  the  trend  of  future  cycles*  The  methods  of  prediction  previously 
used  have  covered  a  wide  range,  including  harmonic  analysis,  assumption 
that  cycles  repeat  after  a  number  of  years,  and  numerous  empirical  rela¬ 
tionships  between  heights  of  maxima,  rate  of  rise,  etc*  Until  recently 
these  predictions  were  largely  of  academic  interest,  but  with  the  dis¬ 
covery  of  the  close  relationships  which  exist  between  radio  propagation 
conditions  and  sunspot  activity  the  prediction  of  sunspot  numbers  has 
assumed  great  practical  importance* 

It  would  be  desirable  to  predict  the  ionospheric  conditions 
directly,  but  that  is  not  yet  possible  since  ionospheric  measurements 
do  not  extend  much  beyond  one  cycle  of  solar  activity — not  enough  to  be 
regarded  as  q.  fair  sample  of  cyclic  phenomena  in  the  statistical  sense. 
However,  the  Z&rich  relative  sunspot  numbers  constitute  a  highly  homo¬ 
geneous  series  extending  back  to  1849,  while  less  reliable  date,  extend 
the  series  back  to  1749*  Since  the  correlation  is  very  high  between  iono*= 
spheric  phenomena,  for  the  interval  during  which  they  have  been  observed, 
and  relative  sunspot  numbers  for  the  same  interval,  conclusions  reached 
from  study  of  the  sunspot  data  may  be  assumed  to  apply  to  the  ionospheric 
phenomena* 

The  various  formulas  which  have  been  developed  for  prediction  of 
the  sunspot  cycle  lead  to  widely  different  predictions.  The  originators 
of  the  various  prediction  methods  do  not  give  satisfactory  estimates  of 
the  reliability  of  their  formulas  for  the  future* 

2.  Method 

To  develop  a  satisfactory  sunspot  prediction  formula  for  use  in 
calculating  rauio  frequencies  for  future  use  the  following  assumptions 
were  made  at  this  laboratory,  following  a  method  developed  by  ttianer: 

(1)  In  a  time  series  exhibiting  cyclic  tendencies  an  estimate,  to  a 
first  approxima tion,  of  a  future  value  in  the  series  is  the  mean  of  all 
past  values  for  the  same  part  of  the  cycle,  and  (2)  this  estimate  can 
be  improved  by  adding  to  the  mean  a  correction  proportional  to  the  de¬ 
partures  of  earlier  values  of  the  same  cycle  from  their  respective 
means,  the  factors  of  proportionality  being  determined  by  the  method  of 
least  squares*  The  prediction  formula  then  becomes 

Hk  *  ~<n  ♦  iKn  =  r‘n  ♦  Vl  Aan-1  *  V-2  a“n-2  ™  U) 

in  which  is  the  value  to  be  predicted  in  a  particular  cycle,  Rn  the 
mean  of  all  the  nth  values  in  preceding  cycles,  &Rn<_£  departure  of 
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the  particular  hn=4  from  and  kn=cj_  is  a  prediction  coefficient  to 

he  calculatedo  The  value  of  Ihe  kn<=£*g  are  adjusted  so  that 

N  0 

Z  (  AHn  » 

is  a  minimum  when  the  summation,  is  taken  over  N  cycles  of  past  data® 

Substituting  for  the  value  of  AR^  in  the  above  expression  its  value 
as  given  in  ( l)  $,  differentiating  with  respect  to  each  and  equating  to 

zero,  we  derive  the  array  of  symmetrical  simultaneous  equations  which  must  be 
solved  to  obtain  the  best  value  of  each  kn=ij_o  They  are 
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in  which  Sn_jt  n.i  =  I  (  A  i^.j  AVRn 

An  array  of  equations  is  formed  for  each  year  of  the  cycle  and  as 
many  coefficients  as  are  desired  may  be  calculatedo  By  arbitrarily  set- 
ting  any  or  several  kn„.ps  equal  to  zero,  coefficients  may  be  obtained 
which  apply  for  predictions  from  particular  years  or  any  combination  of 
themo  Thus,  if  it  is  desired  to  predict,  say,  AE  for  the  fourth  year  of 
the  cycle,  the  coefficients  may  be  derived  for  a  prediction  based  on 
A  R^s  for  the  first,  second  and  third  years,  the  third  year  alone,  the 
second  year  alone  or  any  other  combination »  Naturally  the  coefficient 
for  predicting  from  any  given  year  will  be  different  if  other  years  are 
also  used  in  the  prediction » 

Calculation  of  several  kn„j_ss  was  performed  for  each  year  of  the  sun¬ 
spot  cycle  except  the  zero  years  Data  employed  for  the  computations  were 
the  smoothed  12-month  means  of  Zurich  relative  sunspot  numbers  from  1834 
to  1943°  For  reasons  to  be  set  forth  later,  data  prior  to  1834  were  re- 
jectedo  First,  the  beginning  of  each  cycle  to  the  nearest  tenth  of  year 
was  selected,  each  cycle  being  assumed  to  begin  at  a  point  of  minimum 
sunspot  number  obtained  by  passing  a  smooth  curve  through  the  twelve- 
month  running  averageso  The  times  of  minimum  so  selected  for  the  entire 
series  back  to  1755  are  shown  in  Table  1  compared  with  times  given  by- 
Zurich,  selected  for  being  the  lowest  twelve-month  average  value  for  each 
cycle o  The  differences,  where  they  exist,  are  slighto 

Ha1  ing  selected  the  minimum,  the  twelve-month  means  were  then  ordered 
for  the  successive  years  after  minimum  for  all  the  cycles  and  their  ave¬ 
rages  for  the  last  10  cycles  were  obtainedo  Departures  of  the  twelve- 
month  averages  from  corresponding  values  for  the  mean  cycle  were  then 
derived  for  the  first,  second,  third  year,  etc»,  after  the  minimum  for 
each  cyeleu  These  departures  are  shown  in  Table  2,  together  with  the 
departures  for  the  earlier  cycles  from  the  same  mean  cycle o  These  are  the 
values  of  ARft0 
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Values  calculated  for  Sn_.  n_  ^  are  given  in  Table  3»  Since 

(n-l)-i  =  n-(i  +  l)  etc.  the  number  of  values  of  S  required  for  com¬ 
puting  a  given  number  of  k’s  is  relatively  small*  From  the  values  in 
Table  3,  130  k’s  were  computed  on  the  assumptions  that  A  would  be 
calculated  (l)  from  the  immediately  preceding  value  of  Ad,  (2)  from 
several  successive  preceding  values  of  An,  and  (3)  any  combination 
of  preceding  values  of  An  from  which  certain  values  of  Aa  were  neg¬ 
lected*  In  the  case  of  the  fifth  year  of  the  cycle  values  of  kn_i  were 
computed  for  all  values  of  i  up  to  and  including  5.  By  arbitrarily 
setting  certain  f g  equal  to  zero,  alternate  values  of  the  other 

kfl-jMg  were  computed  which  were  anplicable  for  computing  ARn  from 
specific  preceding  values  ol  A  R.  Thus  for  the  fifth  year  of  the  cycle, 
fifteen  prediction  formulas  were  obtained  based  on  prediction  from  the 
preceding  year,  the  two  preceding  ye&xs,  up  to  the  five  preceding  years 
as  well  as  prediction  from  the  fifth,  fourth,  or  any  other  preceding 
vear  alone  or  any  combination  of  these  years.  The  coefficients  are  given 
in  '"able  4. 

3.  Reliability 

-Predictions  v;ere  made  for  each  of  the  years  of  the  ten  cycles  in 
question  using  all  of  the  combinations  of  prediction  formulas  given  in 
Table  4*  Differences  between  observed  and  predicted  values  were  then 
foi*med  to  test  the  reliability  of  the  predic tions .  An  estimate  of 
reliability  for  predictions  beyona  the  observed  cata  requires  that  the 
stancerd  deviation  (<7- =  ,/£(«'  -  r)2/B)  be  multiplied  by  \/(N/ n-m) in 

which  N  is  the  number  of  cases  and  b,  is  Lhe  number  of  parameters  used 
in  the  prediction,  bhen  the  prediction  is  maae  from  the  mean  alone, 
that  is,  by  assigning  to  R^  the  value  obtained  for  the  nth  year  by 
averaging  all  of  the  cycles,  i«i  =  l;  when  the  value  of  R^  is  predicted 

by  correcting  this  value  by  the  departure  for  the  preceding  year,  M  =  2j 
and  if  the  prediction  is  from  the  five  preceding  years  (as  it  may  be 
when  n  =  5),  =  6.  These  estimates  ol  the  values  of  <J~  for  each  year 

of  the  cycle  axe  given  in  Table  3® 

Kfhen  these  values  of  a~  are  examined,  if  is  apparent  that  for 
most  of  the  years  of  the  cycle  the  "best"  prediction  is  obtained  by 
correcting  the  mean  value  by  the  departure  from  mean  of  the  preceding 
year®  However,  there  are  a  number  of  cases  in  which  the  "best"  predic¬ 
tion  is  obtained  by  using  special  combinations  of  years.  One  might 
naively  assume  that  improved  predictions  would  result  by  selecting  for 
each  year  of  the  cycle  that  particular  group  of  coefficients  which  gave 
rise  to  the  lowest  value  of  <j~  in  the  past  cycles.  This  would  be  falla¬ 
cious,  because,  owing  to  statistical  considers  tions ,  it  is  likely  that 
some  particular  system  may  "work  best"  for  a  limited  number  of  cases. 

The  system  which  "works  best"  on  the  average  should  be  adopted. 

From  the  foregoing  considexa tions  it  was  decided  to  employ  only 
the  mean,  corrected  by  the  departure  of  the  preceding  year,  lor  px-edic- 
tion  of  future  sunspot  numbers.  For  most  practical  considerations  this 
is  satisfactory  for  forecasts  of  propagation  conai tions.  If  more  ad¬ 
vanced  predictions  are  desired,  the  prediction  may  be  based  on  the  mean, 
corrected  for  the  departure  two  years  preceding  the  year  to  be  predicted. 
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Since  the  prediction  coefficients  are  subject  to  statistical 
fluctuations,  it  was  decided  to  plot  them  against  the  year  of  the  cyele 
for  which  they  apply  and  to  adopt  prediction  coefficients  from  a  smoothed 
curve  through  the  values .  By  interpolation  along  this  smoothed  curve, 
prediction  coefficients  may  be  obtained  for  non-integral  times  after  the 
beginning  of  a  cycle,  Fig.  1. 


4»  auto-correlation  in  the  Da ta 


The  calculations  performed  for  obtaining  the  prediction  coefficients 
may  be  employed  for  studying  statistical  characteristics  of  sunspot  cycle. 
A  somewhat  idealized  representation  of  the  data  is  shown  in  Fig.  2k 
obtained  by -spacing  the  minima  shown  in  Table  1  at  11-year  intervals 
and  plotting  the  observed  annual  value  at  these  minima  and  for  the  10 
succeeding  years.  Deviations  of  these  observed  annual  values  from  the 
Mmean  cycle"  obtained  from  the  data  of  the  last  10  cycles  are  shown  in 
Fig.  2B.  The  conservatism  of  the  series  is  conspicuous  from  the  tendency 
of  adjacent  annual  deviations  to  have  the  same  sign  and  similar  magni¬ 
tudes,  the  tendency  upon  which  this  system  of  prediction  is  based 0  It 
is  aleo  clear  that  this  tendency  decreases  as  the  separation  in  time  of 
the  deviations  in  any  given  cyele  is  increased,  so  that  the  further  in 
advance  a  prediction  is  made  the  greater  is  its  uncertainty.  Prediction 
of  the  later  years  of  a  cycle  from  the  early  years  seems  unwarranted. 

These  tendencies  find  numerical  expression  in  the  mean  auto-correla¬ 
tion  coefficients  of  the  series  which  may  be  calculated  from  the  S«v&lues 
by  the  formula 


(j  =  i) 


These  values  are  plotted  in  Fig.  3  for  values  of  i  up  to  3.  It  is 
interesting  to  note  that  if  the  correlation  coefficients  between  each 
pair  of  adjacent  years  in  the  cycle  are  multiplied  by  the  ratios  of  the 
sunspot  numbers  of  the  later  year  to  the  earlier  year  in  the  same  pair, 
a  series  of  prediction  coefficients  is  obtained  which  agree  closely 
with  those  computed  above  from  the  AH’s  alone. 

The  foregoing  remarks  apply  in  particular  to  the  series  in  Fig.  2 
subsequent  to  1834  which  were  used  in  the  calculations.  It  may  be  noted 
that  the  deviations  as  shown  in  Fig.  2B  prior  to  1834  are  somewhat  greater 
than  those  subsequent  to  1834»  There  is  good  reason  for  believing  that 
the  earlier  observations  are  subject  to  considerable  uncertainty  because 
of  absence  of  a  uniform  system  of  observations  and  reductions  from  one 
observer  to  another.  Also  the  data  subsequent  to  1834  are  substantiated 
by  regular  measurements  of  geomagnetic  fluctuations,  initiated  by  Gauss 
and  others  at  that  time,  which  show  a  high  correlation  with  sunspots. 


5 


5e  Test  for  Validity  of'  Data 

It  is  of  interest  to  inquire  whether  any  statistical  basis  for 
rejecting  the  earlier  observations  can  be  established.  In  order  to 
increase  the  size  of  the  working  sample  all  11  years  of  each  of  the 
later  10  cycles  were  normalized  by  expressing  the  deviations  from  the 
mean  cycle  in  terms  of  the  standard  deviation  for  that  year,  thus  sup¬ 
plying  a  sample  of  11  x  10  =  110  individuals.  Frequencies  of  occur¬ 
rence  of  these  deviations  in  6  class  intervals  were  tested  by  the  chi- 
squared  test,  to  ascertain  the  probability  that  they  arose  from  a  normal 
distribution  (degrees  of  freedom  6-3  =  3)«  Neglecting  the  effect  of 
auto-correlation,  there  was  only  1  chance  in  1000  that  the  normalized 
deviations  came  from  a  normal  distribution.  (The  distribution  could 
not  be  expected  to  be  normal  since  the  sunspot  numbers  have  a  lower 
bound  of  zero.) 

The  deviations  were  then  fitted  by  inspection  to  a  type  VI 
Pearsonian  curve  (see  Fig.  4) »  and  the  chi-squared  test  (chi-squared 
=  3°37)  showed  a  probability  of  0.19  that  the  hypothesis  of  a  VI 
Pearsonian  distribution  was  valid,  effects  of  auto- correlation  again 
being  neglected  (in  this  case  the  degrees  of  freedom  were  6-4=2). 

The  observations  prior  to  1834  were  normalized,  using  the  cr  from  the 
later  series,  and  tested  to  ascertain  the  probability  that  they  came 
from  the  same  set  as  the  later  observations.  Chi-squared  in  this  case 
was  117.  In  this  case  the  degrees  of  freedom  were  6-1  =  5,  the  mean, 
standard  deviation,  and  skewness  having  been  determined  by  the  data 
from  the  later  series.  Values  of  chi-squared  as  great  as  117  for  5 
degrees  of  freedom  do  not  appear  in  the  tables;  there  is  practically 
zero  probability  that  the  deviation  obtained  for  the  earlier  years 
of  the  series  came  from  the  same  set  as  those  of  the  later  years. 

The  above  calculations  were  performed  on  the  assumption  that 
there  are  110  independent  observations  in  the  later  series  and  77 
independent  observations  in  the  earlier  series.  A  difficulty  arises 
in  taking  into  account  the  equivalent  number  of  independent  observa¬ 
tions;  certainly  the  numbers  cannot  be  as  great  as  110  and  77,  nor  can 
they  be  less  than  10  and  7.  Fortunately  it  is  not  necessaiy  to  decide 
on  the  exact  equivalent  number  of  independent  observations.  Assuming 
the  extreme  case  that  there  were  10  and  7  independent  observations, 
the  probability  that  the  later  series  comes  from  a  VI  Pearsonian  dis¬ 
tribution  becomes  about  0.35  while  the  probability  that  the  earlier 
series  cemes  from  the  same  distribution  lies  between  0.01  and  0.001. 

This  is  equivalent  to  saying  that  if  we  observe  sunspot  numbers  for  the 
next  10,000  years,  we  would  not  expect  to  encounter  a  series  of  7  cy¬ 
cles  which  differed  as  much  from  the  last  10  cycles  as  the  7  preceding 
them  did.  Since  the  first  auto- correlation  coefficient  is  0.73,  we 
may  assume  that  about  %  of  each  observation  is  independent,  and  that 
therefore  the  number  of  independent  observations  in  the  later  and 
earlier  series  is  110/4  and  77/4  respectively.  These  give  values  of 
chi-squared  equal  to  0.34  and  29  which  indicate  probabilities  of  about 
0.65  and  less  than  0.001  for  the  two  hypotheses. 
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For  the  reasons  stated  above  the  earlier  data  were  not  employed 
in  deriving  the  prediction  formulas®  However,  if  they  had  been  included, 
no  outstanding  change  in  the  coefficients  would  have  resulted  but  the 
estimate  of  reliability  of  the  predictions  would  have  been  deereasedo 

Methods 


One  of  the  faults  of  many  of  the  methods  for  predicting  th® 
sunspot  cycle  is  that  no  reliability  estimate  of  the  accuracy  of  the 
predictions  is  supplied,  a  criticism  which,  it  is  felt*,  cannot  be  made 
against  the  method  set  forth  here®  Estimation  of  reliability  for  csr° 
tain  types  of  formulas  is  very  difficult®  Undoubtedly  many  formulas 
can  be  derived  which,  when  applied  to  past  cycles,  give  satisfactory 
predictions,  for  it  is  always  possible  with  a  sufficiently  large  number 
of  degrees  of  freedom  to  fit  satisfactorily  any  finite  group  of  events 
to  any  desired  accuracy®  In  estimating  the  reliability  of  future  pre¬ 
dictions  due  allowance  must  be  made  for  the  degrees  of  freedom  involved 
in  the  prediction  formula® 

A  more  subtile  consideration  is  involved  in  estimating  the 
reliability  of  a  prediction  based  on  an  ^empirical48  formula®  Empiric 
cal  formulas  are  usually  obtained  by  trying  various  relationships  and 
selecting  those  which  wworkedn  best  in  the  past®  If  a  large  number  of 
relationships  were  tried  it  would  not  be  surprising  to  find  that  some 
of  them  worked  very  well  but  allowance  must  be  mad©  for  the  number  of 
relationships  examined®  Th©  number  implicitly  examined  but  rejected 
by  Simple  inspection  of  the  data  is  ordinarily  much  greater  than  the 
number  explicitly  examined  and  consequently  is  difficult  to  take  into 
account  in  assessing  reliability®  Such  considerations  do  not  apply  if 
an  &  priori  prediction  formula  is  assumed,  even  though  the  parameters 
of  the  formula,  are  obtained  empirically  from  past  data®  But  in  e$ti®&- 
ting  the  reliability  of  the  prediction,  allowance  must  be  made  for  the 
number  of  parameters  in  the  formula  as  was  done  in  computing  the  stan~ 
dard  deviations  of  the  predictions  given  in  Table  5® 

•  7®  Application  to  Present  Cycle 

These  prediction  methods  have  been  applied  to  the  first  part 
of  cycle  18  as  shown  in  Fig®  5  in  which  the  mean  curve  for  the  last 
10  cycles,  the  observed  curve  to  date,  and  the  predicted  points  for  each 
quarter  year  are  represented®  The  predictions  were  made  for  one  year 
in  advarice  using  only  one  coefficient,  naaely  that  to  correct  for  the 
preceding  year®  These  were  interpolated  from  Fig®  1  for  the  quarter^ 
years®  Thus,  the  prediction  for  one  year  after  the  minimum  was  obtained 
by  subtracting  from  the  observed  aainiaura  the  value  of  the  saimimum  for 
the  average  cycle  (7®5  -  5®0  =  2® 5)®  This  quantity  was  multiplied  by 
the  coefficient  relating  the  zero  year  to  the  first  year  (2® 5  S  1®2 
=  3®C,  see  Fig®  1)  and  added  to  th®  value  for  the  first  year  in  the 
average  cycle  giving  the  predicted  value  (3®0  +  15®0  =  18®C)  etc® 
Agreement  wi  th  th©  observed  curve  was  extremely  satisfactory  for  th« 
first  six  quarterly  predictions®  but  rapid  increases  in  the  sunspot 
numbers  during  the  last  few  months  have  thrown  the  predictions  into 
disagreement  for  the  last  available  quarter-year ®  However,  the  dif~ 
ference  between  observed  and  predicted  values  is  only  17® 5  units 
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while  the  standard  deviation  of  the  predictions  is  11®S  units*  There 
is  one  chance  in  six  (approximately)  that  a  value  will  exceed  the 
standard  deviation  by  this  amount® 

The  highest  point  predicted  by  this  method  has  the  value  121.. 

The  limits  of  its  reliability  may  be  stated  as  follows:  Probability 
that  difference  (X)  of  the  observed  and  predicted  values  lies  within 
the  ranges  i>  i  20  units  ,  Go  08;  4-  20  >  X  >  +  10,  0*12;  +  10>X>0, 
Qo23;  0  >  X  >  -  10,  0*32;  -  10  >  X>  -  20,  0*23;  X  >  -  20,  0*03o  Two 
points  are  plotted  on  the  curve  for  more  than  a  year  in  advance;  these 
indicate  that  the  maximum  should  occur  during  the  latter  half  of  1947* 

8®  Comparison  with  Stewart’s  Method 

Plotted  in  Fig*  5  are  also  shown  a  series  of  predictions 
(interpolated  in  part)  by  Prof*  John  Stewart,  communicated  to  this 
laboratory  in  a  note  dated  June  14?  1946*  Prof*  Stewart's  explanation 
of  the  method  of  prediction  is  as  follows:  "A  formula  for  representing 
smoothed  spot  numbers  during  a  single  cycle  has  been  described  in  pre¬ 
vious  papers  (J*  Q*  Stewart  and  H*  A*  Panofsky,  Ap®  J*  88,  pp*  385-407, 
1938;  Jo  Qo  Stewart  and  F*  C*  Eggleston,  Ap*  J*  91,  pp®  72-82,  1940*)* 
When  the  last  previous  cycle,  number  17  had  just  reached  maximum,  at 
the  end  of  1937,  Stewart  and  Eggleston  (J*  Q*  Stewart  and  F®  C*  Eggle¬ 
ston,  Physical  Review  55,  1102,  1939;  and  elsewhere®),  fitted  a  4- 
parameter  curve  to  the  smoothed  monthly  numbers  and  made  a  prediction 
of  the  remaining  course  of  that  cycle,  1938-1944°  This  prediction 
was  successful  within  the  announced  tolerance® 

“Besides  describing  the  4°parameter  representation  (parameters 
s®  a®  fog  F),  Stewart  and  Panofsky  gave  correlations  (p«402,  equations 
(2l!8  {22)  of  b  and  F,  respectively  with  a,  permitting  reduction  to  2 
parameters,,  namely  a  and  a when  the  data  do  not  warrant  more  careful 
fitting®  At  the  time  of  writing,  monthly  numbers  for  the  present  cycle 
are  available  only  through  October,  1946®  The  maximum  is  still  a  year 
ahead,  and  the  limited  part  of  the  curve  yet  observed  is  not  enough  to 
establish  a  really  reliable  prediction*  Nevertheless,  it  is  of  interest 
to  make  the  best  on®  we  can  at  so  early  a  stage  in  the  cycle*  The  solu¬ 
tion  will  improve  with  every  added  month* 

wM.r®  Georg®  S®  Baldwin,  Jr®,  of  Brookline,  Massachusetts,  a 
Princeton  undergraduate  (class  of  1949)?  has  computed  and  plotted  the 
family  sunspot  curves  for  different  values  of  a,  which  follow  when  the 
corresponding  values  of  b  and  F  given  by  the  correlations  mentioned 
above  are  adopted*  To  the  same  scale  monthly  numbers  since  1943  were 
plotted  on  tracing  paper®  By  sliding  the  observed  curve  in  the  time 
coordinate,  a  preliminary  value  of  a  for  a  fit  was  found,  and  also  of  s 
the  time  of  start  of  cycle  18®  Then  the  fit  was  improved  by  a  more 
detailed  comparison  with  carefully  smoothed  monthly  numbers® 

"The  result  — -  which  it  must  b©  emphasized  again  is  only  a 
tentative  fit  for  cycle  18  is  s  =  1943°9?  time  of  start; 


a  —  4*2; 

j  =  1943 .0,  time  of  maximum; 

V  =  233 >  height  of  maximum 

The  computed  curve  reached  spot  number  50  at  1945 * 3 ,  and  is  due  to  reach 
100  at  19466?  and  125  at  1947*3*  The  time  o  maximum  is  subject  to  less 
uncertainty  than  its  height  V  *.  The  predicted  133  would  mean  the  highest 
sunspot  activity  since  1870*  If  spots  during  1946  increase  faster  than 
the  prediction,  the  maximum  will  be  higher  and  will  come  a  tenth-year  or 
more  earlier;  if  the  spots  increase  more  slowly,  the  maximum  will  be  a 
little  later  and  lower* 

"The  values  of  b  and  F  indicated  by  the  correlation  already 
mentioned  and  used  in  the  prediction  ares  b  =  1*027,  log  F  -  lo375° 

"For  convenience  the  formula  for  a  cycle  is  repeated  heres 

ft  =  F(r  -  s)a  e  ~b<r  '  s> 

where  r  is  the  date  in  years  and  R  is  the  computed  (smoothed)  Wolf 
number,  while  e  is  the  base  of  natural  logarithms 0 

"It  should  be  added  that  the  writer  has  by  a  detailed  statistical 
study  established  the  existence  of  an  8-months*  period  in  spot  numbers, 
of  amplitude  sufficient  to  account  for  roughly  half  the  departures  of 
actual  monthly  numbers  from  a  smooth  curve*  (This  period  may  not  be 
rigorously  simple-harmonic,  and  the  phase  may  shift  at  occasional 
intervals  of  a  number  of  yearso)  In  the  present  little  study,  therefore, 
the  observed  monthly  numbers  were  smoothed  by  an  8-month  chain  average*" 

Profo  Stewart’s  note  was  in  hand  when  the  above-described  investi- 
gallon  was  begun*  The  essential  difference  between  the  two  methods  lies 
in  the  fact  that  his  method  makes  no  allowance  for  regression  toward  the 
mean  cycle*  To  date  Prof*  Stewart's  predictions  are  somewhat  better 
than  those  derived  by  the  more  elaborate  method  described  in  this  report 
but  it  may  be  expected  that  as  the  present  sunspot  cycle  develops  the 
law  of  averages  will  operate  in  favor  of  the  statistical  method* 
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Table  1. Years  of  sunspot  minims  for  cycles  1-18 


Cycle  no. 

Selected  for 
this  study 

Selected  by 
Zurich 

1 

1755.2 

1755.2 

•> 

1766.5 

17b6. 5 

3 

1775.6 

1775.5 

4 

1734.6 

1784.7 

- 

l'K!3.3 

1798.3 

A 

1  S3. 0.6 

1310.6 

7 

1823.2 

1823.3 

8 

1834.1 

1833.9 

9 

1843.6 

1843.5 

10 

1856.3 

1856.0 

13 

1867.3 

1867.2 

12 

1578.8 

1878.9 

13 

1889.5 

1389.6 

H 

1903.8 

1901.7 

15 

1933.4 

1 913 „ ' 

16 

1923.6 

1923.6 

17 

1033.3 

1933.8 

18 

1944.2 

1°44.  ' 

Table  3. Cross  product  summsticn  terns  {Sn„-\ , 
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coefficient  set  equal  to  zero  in  computation. 


Table  5. Standard  deviations  of  predictions  for  cycles  8-17  using  coefficients,  kn.^,  of  Table 


b 


g  I.  VALUE  OF  PREDICTION  COEFFICIENTS,  k,  USING  ONLY  PRECEDING  YEAR’S 
DATA  FOR  PREDICTING  DEPARTURE  (  A  Rn  =  K  A  Rn  _  ,  ) . 


Fig.  2A.  ELEVEN  YEAR  SERIES  OF  ANNUAL  SUNSPOT  Nl  ,V,BERS  1834-1944  (A)  AND  DEVIATION  FROM  MEAN  CYCLE  (B) 
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Fig.  3.  CORRELOGRAM  SHOWING  AUTO- CORRELATION  COEFFICIENTS  (r 
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Fig.  4.  PEARSON  YL  CURVE  FIT  TO  YEARLY  DEVIATIONS  FROM  MEAN  CURVE  FOR  CYCLES  8-17 
AND  CYCLES  I -7. 
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Fig.  5.  COMPARISON  OF  OBSERVATIONS  AND  PREDICTIONS  FOR  CYCLE  lb. 


